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Abstract 



We present results on the production of high transverse momentum tt and 77 
mesons in pp and pBe interactions at 530 and 800 GeV/c. The data span the 
kinematic ranges: 1 < pT < 10 GeV/c in transverse momentum and 1.5 units 
in rapidity. The inclusive 7r° cross sections are compared with next-to-leading 
order QCD calculations and to expectations based on a phenomenological 
parton-A;^ model. 

PACS numbers: 12.38.Qk, 13.85.Ni 
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Typeset using REVT^ 



I. INTRODUCTION 



The study of inclusive single-hadron production at large transverse momentum (px) 
has been a useful probe in the development of perturbative quantum chromodynamics 
(PQCD) [1,2]. Early in the evolution of the parton model, a departure from an expo- 
nential dependence of particle production at lower px was interpreted in terms of the onset 
of interactions between pointlike constituents (partons) contained in hadrons. Large px is 
a regime where perturbative methods have been applied to QCD to provide quantitative 
comparisons with data. Such comparisons yield information on the validity of the PQCD 
description, and on parton distribution functions of hadrons and fragmentation functions of 
partons. 

This paper reports high-precision measurements of the production of 7r° and rj mesons 
with large px- The 7r° production cross sections are compared with next-to-leading order 
(NLO) PQCD calculations [3]. As illustrated in a previous publication [4], our data (for both 
inclusive tt^ and direct-photon production) are not described satisfactorily by the available 
NLO PQCD calculations, using standard choices of parameters. Similar discrepancies have 
been observed [5] between conventional PQCD calculations and other measurements of high- 
Px 7r° and direct- photon cross sections (see also [6-8]). The origin of these discrepancies may 
be attributed to the effects of initial-state soft-gluon radiation. Such radiation generates 
transverse components of initial-state parton momenta, referred to below as kx [9] • Evidence 
of significant kx in various processes, and a phenomenological model for incorporating its 
effect on calculated high-py cross sections, have been extensively discussed in Refs. [4,5]; 
recent studies of the photoproduction of direct photons at HERA may provide additional 
insights [10-14]. The inadequacy of NLO PQCD in describing /cr-sensitive distributions has 
been discussed in Ref. [15]. In this paper, we follow the phenomenological prescription of 
Ref. [5] when comparing calculations with our 7r° data. We also present cross sections for 
T) meson production at large-pT- As might have been expected from earlier measurements 
(see, e.g., [16]), rj production, relative to 7r° production, shows little dependence onp^ or on 
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center of mass rapidity (ycm)- 



II. THE EXPERIMENTAL SETUP 

The E706 experiment at Fermilab was designed to measure direct-photon production 
at high-pr, and to investigate the structure of events containing direct photons. The data 
collection phase of the experiment spanned three fixed-target running periods, and included 
a relatively low statistics commissioning run in 1987-88 [16-19], and primary data runs in 
1990 and 1991-92. The results presented here are from data recorded during the 1991-92 run. 
The E706 apparatus, operated in tandem with the E672 muon spectrometer, constituted the 
Meson West Spectrometer, displayed schematically in Fig. 1. The experiment used a right- 
handed Cartesian coordinate system, with the Z-axis pointed in the nominal beam direction, 
and the Y-axis pointed upward. The principal elements of the Meson West Spectrometer 
are discussed below. A more detailed description of the triggering methodology appears in 
Ref. [20]. 

A. Beamline and target 

The Meson West beamline was capable of transporting either 800 GeV/c primary protons 
from the Fermilab Tevatron or secondary beams of either polarity. We report here on results 
from studies using 800 GeV/c primary protons and a 530 GeV/c positive secondary beam. 
The beamline was instrumented with a differential Cerenkov counter [21,22] to identify 
incident pions, kaons, and protons in secondary beams. This helium-filled counter was 
43.4 m long and was located ^100 m upstream of the experiment's target. Using the 
Cerenkov counter, the proton fraction at 530 GeV/c was determined to be 97% [21]. 

A 4.7 m long stack of steel surrounding the beam pipe was placed between the last 
beamhne magnet and the target box (see Fig. 1) to absorb hadrons. A water tank was placed 
at the downstream end of this hadron shield to absorb low-energy neutrons. During the 1991- 
92 run, two walls of scintillation counters were located both upstream and downstream of 
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the hadron shield to identify penetrating muons. 

The target region during the 1991-92 run consisted of two 0.8 mm thick copper disks of 
2.5 cm diameter, located immediately upstream of a liquid hydrogen target, followed by a 
2.54 cm long beryUium cylinder of diameter 2.54 cm. The hydrogen target consisted of a 
15 cm long mylar flask, supported in an evacuated volume, with beryllium windows at each 
end (2.5 mm thickness upstream and 2.8 mm thickness downstream) [23]. 

B. Charged particle tracking 

The spectrometer employed a charged particle tracking system consisting of silicon strip 
detectors (SSDs) [24], a dipole analysis magnet, proportional wire chambers (PWCs), and 
straw tube drift chambers (STDCs) [25]. The SSD system consisted of sixteen planes of 
silicon wafers, arranged in eight modules. Each module contained two SSD planes, one 
providing X-view information and the other instrumenting the F-view. Six 3x3 cm^ SSD 
planes were located upstream of the target and used to reconstruct beam tracks. Two hybrid 
5x5 cm^ SSD planes (25 fj,m pitch strips in the central 1 cm, 50 beyond) were located 
Ki2 cm downstream of the Be target [26]. These were followed by eight 5x5 cm^ SSD planes 
of 50 iim pitch. The SSDs were instrumented to cover a solid angle of ±125 mr. Figure 2 
displays the distribution of reconstructed vertices from a representative sample of 1991-92 
data, showing clear separation of the different target elements. 

The analysis dipole magnet imparted a transverse momentum impulse of Ri450 MeV/c (in 
the horizontal plane) to singly-charged particles. Downstream track segments were measured 
by means of four stations of four views (XYUV) of 2.54 mm pitch PWCs and two stations of 
eight (AXAY) layers of STDCs with tube diameters 1.03 cm (upstream station) and 1.59 cm 
(downstream station). The STDC stations were installed prior to the 1990 flxed-target run, 
and improved the angular resolution of the downstream tracking system (to paO.OG mrad) to 
make it comparable to that of the upstream system. 
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C. Calorimetry 



The central element of the E706 apparatus was the finely-segmented hquid argon electro- 
magnetic calorimeter (EMLAC) used to detect and measure electromagnetic showers [27,28]. 
The EMLAC had a cylindrical geometry with an inner radius of 20 cm and an outer radius 
of 160 cm. It was divided into four mechanically independent quadrants, which were further 
subdivided electronically to create octants. The calorimeter had 33 longitudinal cells read 
out in two sections: an 11 cell front section {^8.5 radiation lengths) and a 22 cell back 
section (?»18 radiation lengths). This front/back split was used for measuring the direction 
of incidence of showering particles, for discriminating between electromagnetic and hadronic 
showers and for resolving closely separated electromagnetic showers. The longitudinal cells 
consisted of 2 mm thick lead cathodes (the first cathode was constructed of aluminum), 
double-sided coppcr-clad G-10 radial (R) anode boards, followed by 2 mm thick lead cath- 
odes and double-sided copper-clad G-10 azimuthal ($) anode boards. There were 2.5 mm 
argon gaps between each of these layers in a cell. The physical layout is illustrated in Fig. 3. 

The copper-cladding on the anode boards was cut to form strips. Signals from corre- 
sponding strips from all R (or $) anode boards in the front (or back) section were jumpered 
together. The copper-cladding on the radial anode boards was cut into concentric strips 
centered on the nominal beam axis. The width of the strips on the first R board was 
5.5 mm. The width of the R strips on the following R boards increased slightly so that the 
radial geometry was projective relative to the target, which was located 9 m upstream of the 
EMLAC. The azimuthal readout was subdivided at a radius of 40 cm into inner and outer 
segments, with each inner $ strip subtending an azimuthal angle of 7r/192 radians, and outer 
$ strips covering 7r/384 radians. Subdivision of the azimuthal strips in the outer portion of 
the detector improved both the position and energy resolution for showers reconstructed in 
this region. It also reduced correlation ambiguities from multiple showers in the same 
octant of the calorimeter. 

Data acquisition and trigger-signal processing for the EMLAC was based upon the FNAL 
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RABBIT system [29] . To achieve tolerable deadtimes, the zero suppression features of this 
system were used extensively during the experiment's 1987-88 commissioning run. How- 
ever, zero suppression limited our ability to understand the effects of out-of-time photon- 
induced showers and the tails of hadron-induced showers, thereby compromising efforts to 
understand the detailed response of the detector. Consequently, FASTBUS modules (the 
ICBM — Intelligent Control and Buffering Module — and the Wolf interface [30,31]) were de- 
veloped by E706 to replace the original, CDF-designed, MX readout controllers. These 
FASTBUS modules enabled us to eliminate zero suppression during the experiment's two 
primary data runs. 

The apparatus also included two other calorimeters: a hadronic calorimeter (HALAC) lo- 
cated downstream of the EMLAC in the same cryostat, and a steel and scintillator calorime- 
ter (FCAL), positioned further downstream, to increase coverage in the very forward region. 
The HALAC had 53 longitudinal cells read out in 2 sections: a 14 cell front section (?»2 
interaction lengths), and a 39 cell back section (Ri6 interaction lengths). Each cell consisted 
of read-out planes separated by 3 mm argon gaps and a 2.5 cm thick steel plate. 

The FCAL acceptance covered the beam hole region of the EMLAC. It was split into 
three longitudinally similar sections. Each section was composed of alternating layers of 
1.9 cm thick steel absorber plates and 4.8 mm thick acrylic scintillator sheets. The distance 
between steel plates was 6.9 mm. The downstream module contained 32 steel absorber 
plates and 33 scintillator sheets; the other two modules were comprised of 28 steel absorber 
plates and 29 scintillator sheets. Together, the three modules constituted KilO.5 interaction 
lengths of material. 

D. Muon identification 

The E672 muon spectrometer, consisting of a toroidal magnet, shielding, scintillators, 
and proportional wire chambers, was deployed immediately downstream of the FCAL. The 
combined Meson West Spectrometer was triggered on high-mass muon pairs in order to 
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investigate the hadroproduction of J/'^, ip{2s), Xc and B mesons [32-35]. E706 and E672 
collected data simultaneously and shared trigger logic, the data acquisition system, and 
event reconstruction programs. Data collected with the dimuon trigger were also used for 
several technical studies — for example, the J/ip — > A* '''A*" signal was used to calibrate the 
momentum scale of the tracking system. 

E. Triggering 

The E706 trigger selected events yielding high transverse momentum showers in the EM- 
LAC. The selection process involved four stages: beam and interaction definitions, pretrigger 
requirements, and the final trigger requirements. Beam particles were detected using a ho- 
doscope consisting of three planes (arranged in X, Y and U views) of scintillator, located 
Ri2 m upstream of the target region. Each plane contained 12 scintillator strips. The widths 
of the strips varied from 1 mm in the central region, to 5 mm along the edges. The edges of 
individual strips in each plane were overlapped to avoid gaps in acceptance. A BEAM signal 
was generated by a coincidence of signals from counters in at least two of the three hodoscope 
planes. A beamI signal required that less than two hodoscope planes detected two or more 
isolated clusters of hits in coincidence with BEAM. This beamI requirement rejected events 
with multiple beam particles incident upon the target. A plane of four scintillation counters 
(referred to as beam hole counters) arranged to produce a 0.95 cm diameter central hole 
was located downstream of the beam hodoscope and used to reject interactions initiated by 
particles in the beam halo. 

Two pairs of scintillation counters were mounted on the dipole analysis magnet, one 
pair upstream and the other downstream of the magnet. Each pair had a central hole that 
allowed non-interacting beam particles to pass through undetected. An interaction was 
defined as a coincidence between signals from at least two of these four interaction counters. 
To minimize potential confusion in the EMLAC due to out-of-time interactions, a filter was 
used to reject interactions that occurred within 60 ns of one other. 
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For those interactions that satisfied both the beamI and interaction definitions, the 
Pt deposited in various regions of the EMLAC was evaluated by weighting the energy signals 
from the fast outputs of the EMLAC R channel amplifiers by aisin6'i, where 9i is the polar 
angle that the i*'* strip subtends with respect to the nominal beam axis. The pretrigger hi 
requirement was satisfied when thepT detected either in the inner 128 i? channels or the outer 
R channels of any octant was greater than the threshold value of ~2 GeV/c. A pretrigger 
signal was issued only if the signals from a given octant satisfied the pretrigger requirement, 
there was no evidence in that octant of substantial noise or significant pr attributable to 
an earlier interaction, and there was no incident beam halo muon detected by the walls 
of scintillation counters surrounding the hadron shield upstream of the spectrometer. The 
pretrigger signal latched data from the various subsystems until a final trigger decision was 
made. 

Localized trigger groups were formed for each octant by clustering the R channels into 
32 groups of 8 channels. Each of the adjacent pairs of 8 channel groups (groups 1 and 2, 2 
and 3, . . . , 31 and 32) formed a local group of 16 strips. If the pt detected in any of these 
groups of 16 was above a specified high (or low) threshold, then a LOCAL Hi (or LOCAL LO) 
signal was generated for that octant. The SINGLE LOCAL Hi (single LOCAL LO) trigger 
required a LOCAL Hi (local lo) signal from an octant that also satisfied the pretrigger 
HI. The LOCAL HI (local lo) threshold was ^3.5 GeV/c {^2 GeV/c). The SINGLE LOCAL 
LO trigger was prescaled by a factor of Ri200. 

In addition to these high-p^^ triggers, prescaled samples of low-bias beam, interaction, 
and pretrigger events were also recorded. The prescale factors for these triggers were typically 
set at 15^, 15^, and 15^, respectively. These low-bias triggers constituted ^10% of the events 
recorded. 
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III. 7r° AND 7/ ANALYSIS 



The data sample used in this analysis corresponds to an integrated luminosity of 6.8 
(1.2) cvents/pb for 530 GeV/cpBe (pp) interactions, and 6.5 (1.2) events/pb at 800 GeV/c. 
The following subsections describe the analysis procedures and methods used to correct the 
data for losses resulting from inefficiencies and selection biases. 

A. Event reconstruction 

The charged-track reconstruction algorithm produced track segments upstream of the 
dipole magnet using information from the SSDs, and downstream of the magnet using in- 
formation from the PWCs and STDCs. These track segments were projected to the center 
of the magnet, and linked to form the final tracks, and the interaction vertex. The charged- 
track reconstruction and vertex-finding methodology are described in Ref . [20] . 

The readout in each EMLAC quadrant consisted of four regions: left R and right R 
(radial strips of each octant in that quadrant), and inner $ and outer $ regions (azimuthal 
strips divided at = 40 cm). Strip energies from clusters in each region were fitted to 
the shape of electromagnetic showers, as determined from Monte Carlo simulations and 
isolated-shower data. These fits were used to evaluate the positions and energies {Er and 
E^) of the peaks in each region. Shower positions and energies were obtained by correlating 
peaks of approximately the same energy in the R and $ regions within the same half octant 
(more complex algorithms were used to handle configurations with overlapping showers in 
either the it! or $ regions). Any differences in photon energy as measured in the R and $ 
views reflect the intrinsic resolution properties of the calorimeter, and provide a test of the 
quality of the Monte Carlo simulations. The EMLAC's longitudinal segmentation provided 
discrimination between showers generated by electromagnetically or hadronically interacting 
particles. For individual showers, the ratio of energy reconstructed in the front section to the 
sum of energy in the front and back section of the EMLAC (referred to as £^front/-E'TOtal) 
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also tested the Monte Carlo simulation of longitudinal shower development (see the detector 
simulation section below) . An expanded discussion of the EMLAC reconstruction procedures 
and performance can be found in Ref. [28]. 

B. Data sample selection and corrections 

The events contributing to the cross section measurements were required to have a recon- 
structed vertex within the fiducial volume of the Be or H2 targets [36] . Vertex reconstruction 
efficiencies were evaluated for each target using a detailed Monte Carlo simulation of the 
spectrometer (described below). These efficiencies were used to correct for reconstruction 
losses and for resolution smearing across target fiducial boundaries. The vertex reconstruc- 
tion efficiency was 1 for the H2 and downstream Be targets and 0.97 for the upstream Be 
target. 

Both 7r° and r] mesons were reconstructed via their 77 decay modes. The photons were 
required to be within the fiducial region of the EMLAC, which excluded areas with reduced 
sensitivity. In particular, photons incident upon regions of the detector near quadrant 
boundaries (which abutted steel support plates), the central beam hole, the outer radius 
of the EMLAC, and octant boundaries were excluded from consideration. A simple ray- 
tracing Monte Carlo program was employed to determine the correction for the effect of the 
azimuthal fiducial boundaries [37]. 

To reduce backgrounds from hadronic showers, only showers with at least 20% of their 
shower energy deposited in the front part of EMLAC (-Bfront/-Etotal> 0.2) were accepted 
as photon candidates. To correct the cross sections for this requirement, as well as for other 
larger effects including resolution-smearing and reconstruction losses, a full simulation of the 
showers in the calorimeter (described below) was employed. In addition, 77 combinations 
were considered as 7r° or r) candidates only when the two photons were detected in the same 
octant (to simplify subsequent analysis of the trigger response), and only those combinations 
with energy asymmetry [A^^ = \E^-^ — E^^\/{E^^ + E^^)] less than 0.75 were considered (to 
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reduce uncertainties due to low energy photons). 

The 77 invariant mass distribution in the 7r° and rj regions for photon pairs that satisfied 
the above requirements is shown in Fig. 4 for representative low pt intervals. A 7r° candidate 
was defined as a combination of two photons with invariant mass, M^j [38], in the range 
100 MeV/c^ < My^ < 180 MeV/c^. An 77 candidate was defined as a two- photon combination 
in the range 450 MeV/c^ < My^ < 650 MeV/c^. To evaluate production cross sections, 
combinatorial background in the 7r° and 77 regions was evaluated as follows. Sideband regions 
were defined to cover an equivalent mass range of the 7r° and rj peak regions (using the same 
acceptance criteria as for the peak regions). The pr and rapidity distributions from these 
side bands were then subtracted from the corresponding distributions within the 7r° and 77 
mass ranges to obtain the respective signals. This technique is appropriate as long as the 
combinatorial background depends approximately linearly dependence on M^^. At low pt 
(below 2 GeV/c), the shape of the combinatorial background in the signal regions is not 
linear, and a more sophisticated fitting procedure was used to evaluate the background. 
The 77 mass distributions were fitted to Gaussians for signal plus second-order polynomials 
in My^ to represent the background. The combinatorial background in the peak regions 
was then defined using the resultant parameters of the fit, and the signals defined as the 
differences between the totals and the fitted backgrounds. 

For the cross section measurements, the signals have been corrected for losses due to 
the energy asymmetry cut and the branching fractions [39] for the 77 decay modes. The 
correction for losses due to the conversion of one or both of the photons into e^e~ pairs was 
evaluated by projecting each reconstructed photon from the event vertex to the reconstructed 
position in the EMLAC. The number of radiation lengths of material traversed along the 
photon path was calculated on the basis of a detailed description of the detector. The photon 
conversion probability was then evaluated and used to account for losses from conversion. 
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C. Trigger response 



As mentioned previously, the SINGLE LOCAL Hi and SINGLE LOCAL LO trigger decisions 
were based upon depositions of pt in the EMLAC within groups of 16 contiguous radial 
strips. For each 7r° or rj candidate, a probability to satisfy the trigger was defined, based 
upon energy deposition in the entire octant, as P = 1 — n(l ~Pi), where Pi is the efficiency 
of the i*'* trigger group in the octant containing the candidate [40]. The inverse of this 
probability was applied as a trigger weight to each meson candidate. To avoid excessively 
large trigger weights, meson candidates with trigger probabilities of < 0.1 were excluded 
from further consideration. The correction for losses from this requirement was determined 
from Monte Carlo, and absorbed into the reconstruction efficiency. 

The efficiency of the pretrigger hi was determined in a manner similar to that used 
for the local triggers [40]. Additional details on the trigger can be found in Refs. [20,41,42]. 

The cross sections presented in this paper reflect composite measurements, utilizing a 
combination of results from the interaction, pretrigger hi, single LOCAL LO, and 
SINGLE LOCAL HI triggers. The pr spectra (corrected for only prescale factors) from a 
representative sample of these triggers are shown in Fig. 5. The transition points chosen 
between the high and low threshold triggers were determined by comparing the fully cor- 
rected results from each trigger, and were different for 7r° and rj mesons and also depended 
on rapidity. 

D. Rejection of beam halo muons 

Spurious triggers were produced by muons in the beam halo that radiated energy in the 
electromagnetic calorimeter in random coincidence with an interaction in the target. Partic- 
ularly in the outer regions of the EMLAC, this energy was recorded as a high-p^- deposition 
that satisfied the local trigger requirements. This occurred much more frequently in data 
from the 530 GeV/c secondary beam than for the primary 800 GeV/c beam, due to the 
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absence of an upstream interaction target in the latter case. To reduce this background, 
the pretrigger logic relied on signals from the veto walls of scintillator counters to reject 
events associated with such muons in the beam halo. In the off-line analysis, we employed 
expanded requirements on the latched veto- wall signals, the direction of reconstructed show- 
ers [28], the shower shape (halo muon-induced showers have a different shape than photon 
or electron-induced showers), and the total pr (balance) in the event. For the latter, we 
calculated the net pr of the photons and charged particles which, based upon their initial 
trajectories, would intercept the EMLAC in the transverse plane within the 120° sector 
opposite the meson candidate. In interactions which generate a high-p^ meson, the ratio 
of this "away-side" pr (P^""^^) to that of the meson pt should be near unity. However, for 
events triggered by showers from beam halo muons, P^^'^^ /pr should be near zero, since 
the interaction in random coincidence with the beam halo muon is typically a soft (low pt) 
interaction. Candidates with P^^^^/pt < 0.3 were considered likely to be due to beam halo 
muons and were rejected. 

To illustrate the effect of the above off-line requirements, the 77 invariant mass distri- 
bution, both before and after application of the rejection criteria, is shown in Fig. 6 for the 
530 and 800 GeV/c data, for pt > 7.0 GeV/c. The large muon-induced background at low 
77 mass values in the 530 GeV/c data is due to the occasional transverse splitting of the 
muon-induced showers into two closely-separated photon candidates. This happens because 
the reconstruction software assumes that the showers originate from the target-region rather 
than from the beam halo. The 800 GeV/c data have very few muon-induced triggers, and 
is consequently not affected very strongly by these rejection criteria. 

The impact of these rejection criteria on the physics signal was checked using more 
restrictive selection criteria to define a pure sample of 77 pairs. The fraction of signal lost 
by the application of each of the muon-rejection requirements determined a correction to 
the cross section. The product of the correction factors for muon rejection corresponds to 
an increase of ^8% in the cross secton at = 4 GeV/c, and ~10% at = 7 GeV/c, for 
the 530 GeV/c beam data. (The corrections were smaller for the 800 GeV/c data because 
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of cleaner beam conditions.) 



E. Detector simulation 

The Meson West spectrometer was modeled using a detailed geant [43] simulation. 
Because the full simulation of electromagnetic showers requires extensive computing time, 
we developed a hybrid approach using GEANT-tracking through the magnetic spectrometer 
and in the initial stages of the shower development in the calorimeter. We used the standard 
GEANT algorithms for tracking particles with energies above 10 MeV, below which we rehed 
on an empirical parameterization for the deposition of energy in the EMLAC [44]. This 
cutoff was selected to be at the point at which bremsstrahlung still dominates the energy 
loss in lead, and led to significant improvement in processing speed. In doing this, we took 
advantage of the steady advances in computational power of the FNAL UNIX farms [45] to 
reach the desired level of statistical accuracy. 

As inputs to the GEANT simulation, we employed single particle distributions, recon- 
structed data events, and two physics event generators: HERWIG [46] and pythia [47] . For 
the analyses described herein, we chose HERWiG as the principal Monte Carlo event gener- 
ator based on a better match of particle multiplicities between data and Monte Carlo using 
the default parameters. Over 5.5 million HERWiG events were passed through the GEANT 
simulation. We weighted the HERWIG 7r° and r] spectra (in pt and rapidity) to match our 
measured results, so that the corrections obtained from the Monte Carlo were based on the 
data distributions rather than on the behavior of the physics generator. 

The calibration of the energy response of the EMLAC was based on the reconstructed 
masses of 7r° mesons in the 77 decay mode [28]. The steeply falling pt spectrum for 7r° 
production, combined with the calorimeter's resolution, produced a small offset (~1%) in the 
mean reconstructed photon energies. Using the same calibration procedure in the simulated 
EMLAC as in the detector, we corrected this offset and minimized any potential biases in 
the cahbration. We also employed the GEANT Monte Carlo simulation to evaluate the mean 
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correction (as a function of photon energy) for energy deposited in the material upstream 
of the EMLAC. The impact of detector resolution on the energy scale and on the 7r° and 77 
production spectra was incorporated in the overall reconstruction efficiency corrections. 

To ensure that the Monte Carlo simulation reproduced the data, a special preprocessor 
was used to convert GEANT information into signals and strip energies as measured in the 
various detectors, and to simulate hardware effects, such as channel noise and gain variations. 
The generated Monte Carlo events were then processed through the same reconstruction 
software used for the analysis of data, and thereby provided measures of the inefficiencies 
and biases of the reconstruction algorithms. 

Comparisons between results from the Monte Carlo simulation and the data for the 
distributions in Eji — and in the fraction of energy deposited in the front section of 
the EMLAC, £^front/-£'TOtal, are presented in Figs. 7 and 8. The Monte Carlo results 
are in satisfactory agreement with the data, indicating that the simulation properly treats 
shower development in the EMLAC. Figures 9 and 10 show the 77 mass spectra in the 
7r° and r] mass regions for two minimum-p^ requirements, and compare these spectra to 
the Monte Carlo simulation. In addition to giving further evidence that the Monte Carlo 
provides a good simulation of the resolution of the EMLAC, the agreement in the levels 
of the combinatorial background indicates that the Monte Carlo also provides reasonable 
simulation of the underlying event structure. Figure 11 shows a comparison between the 
Monte Carlo and the data for the sideband-subtracted tt*^ energy asymmetry distribution, 
and the agreement indicates that the Monte Carlo simulation describes accurately the losses 
of low-energy photons. 

Reconstruction inefficiencies for 7r° and 77 mesons (which satisfied the A-y^, energy asym- 
metry, requirement) were relatively small over most of the kinematic range. Figure 12 shows 
the probability for a t:^ to pass the selection requirements imposed on the Monte Carlo 
events at 530 GeV/c, as a function of pt, for different rapidity intervals. This probability 
includes losses due to the reconstruction algorithm, the -EfrontZ-^total requirement, and 
the 10% minimum trigger probability requirement. The inefficiency at forward rapidities 
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and high-pT is attributable to the increased difficulty in separating two photons from 7r° 
decays (coalescence) in this kinematic region. 

F. Normalization 

Electronic scalers that counted signals from the beam hodoscope, interaction counters, 
and beam hole counters were used to determine the number of beam particles incident on 
the target. Other scalers logged the state of the trigger and of components of the data 
acquisition system. Information from these scalers was used to determine the number of 
beam particles that traversed the spectrometer when it was ready to record data. This 
number was corrected for multiple occupancy in the beam hodoscope (beyond that excluded 
via the beamI requirement, and for absorption of beam in the target material. 

The normalization of the low pr 7r° cross section was independently verified using events 
from the prescaled beam and interaction trigger samples. In these samples, the absolute 
normalization can be obtained just by counting events. For these low pt events, the nor- 
malization as determined via the scalers and via event counting techniques agree to better 
than 3%. 

Based upon the good agreement between results from these independent normalization 
methods, combined with the stability of the cross section results from various sections of 
the run, an evaluation of the internal consistency of the scalers, and a detailed analysis of 
the design, implementation and performance of the trigger, the net systematic uncertainty 
in the overall normalization is ~8%. 

G. Secondary Beam Contamination 

The 530 GeV/c secondary beam cross section measurements were corrected for the small 
admixture of pions and kaons present in the beam, estimated as 2.75% tt"^ and 0.5% K'^ , 
respectively [21]. Although the percentage contamination is small, it's effect at high pt 
is enhanced by the presence of two, rather than three, valence quarks in incident mesons. 
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The effect of tt"*" contamination was estimated using our high statistics study of 7r° and rj 
meson production by 515 GeV/c 7r~ beam [4,48]. This is justified because: neutral meson 
production by 7r+ and 7t~ beams is expected to be similar, based on isospin arguments and 
earlier measurements [49]; the tt"*" component of the nominally 530 GeV/c positive secondary 
beam had a mean momentum of 515 GeV/c; and the ratio of the measured high pt 7r° cross 
section from the Cerenkov-tagged tt"*" component of the positive secondary beam to the 
corresponding cross section from the 515 GeV/c 7r~ beam was consistent with unity [21]. 

The effect of K'^ contamination was assumed to be half that of the same amount of 
TT"*" contamination, consistent with earlier, lower energy measurements [50] and with our 
own, more statistically limited, data [21]. After correcting for beam contamination, the 
530 GeV/c cross sections were reduced by 2% at low pt and by 10% at high pt- 

H. Summary of systematic uncertainties 

The principal contributions to the systematic uncertainty arose from the following 
sources: calibration of photon energy response (5-9%), 7r° and rj reconstruction efficiency and 
detector-resolution unsmearing (5%), the overall normalization (8%), and, for the 530 GeV/c 
secondary beam, the beam contamination (0-7%). A more complete list of systematic uncer- 
tainties is presented in Table I. Some of these uncertainties {e.g. normalization) are strongly 
correlated between bins. The systematic uncertainties, combined in quadrature, are quoted 
with cross sections in the appropriate tables. Combining all systematic uncertainties yields 
a net uncertainty that varies from 11% (13%) at lowpT for 7r° (77) mesons to 15% at high-p^. 

The secondary proton beam was determined to have a mean momentum of 530 ± 2 GeV/c 
with an estimated halfwidth of ~ 30 GeV/c. This momentum spread introduces a small 
uncertainty (~ 5%) in comparisons of theory with data. (For the 800 GeV primary beam, 
the momentum bite and the corresponding uncertainty are negligible.) 
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IV. INCLUSIVE CROSS SECTIONS 



A. 7r° production 

The inclusive tt*^ cross sections per nucleon versus pr are shown in Figs. 13 and 14 
for 530 and 800 GcV/c proton beams, respectively, incident upon beryllium and Uquid 
hydrogen targets. Because of the steeply falling spectra, the data are plotted at abscissa 
values that correspond to the average values of the cross section in each pt bin (assuming 
local exponential pt dependence) [51]. These cross sections are also tabulated in Tables II 
and III. The corresponding cross section measurements, as functions of pr and t/cm are 
reported in Tables lY through YII. 

NLO PQCD calculations [3] arc compared to the data in Figs. 15 through 23. In Fig. 15, 
NLO PQCD results using CTEQ4M parton distribution functions [52] and BKK fragmen- 
tation functions [53] are compared to the measured inclusive 7r° cross sections for pBe and 
pp interactions at 800 GeV/c. The PQCD calculations for the Be target have been ad- 
justed to account for nuclear effects by the factor A"~^, where A is the atomic number and 
a — 1.12 (1.08) for 530 (800) GeV/c incident protons. Theoretical results are presented 
for three values of the factorization scale: n = J9t/2, Pt, and 2pT- (In these comparisons, 
the renormalization and fragmentation scales have been set to the value of the factorization 
scale.) In addition to a large dependence on the choice of scale, the expectations for these 
choices of fi lie significantly below the data, at both 530 and 800 CcV/c. 

In Fig. 16, NLO calculations using BKK and KKP [54] fragmentation functions, and /i 
— pr/'i, are compared to the n° cross sections for 530 GeV/c pBe and pp interactions. The 
calculations exhibit considerable dependence on the choice of fragmentation function, but 
both choices predict yields that are significantly lower than the data. 

These discrepancies have been interpreted [4-6] as arising from additional soft-gluon 
emission in the initial state that is not included in the NLO calculation, and which results 
in sizeable parton kr before the hard collision (for a different perspective, see the discus- 
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sion in Ref. [8]). Soft-gluon (or kr) effects are expected in all hard-scattering processes, 
such as the inclusive production of jets, high-p^ mesons, and direct photons [55-58]. The 
CoUins-Soper-Sterman resummation formalism [59] provides a rigorous basis for understand- 
ing these radiation effects, and there have been several recent efforts to derive resummation 
descriptions for the inclusive direct-photon [60-63] and dijet cross sections [64-66]. The 
calculation of Ref. [60] for inclusive direct-photon production, which includes the effects 
of soft-gluon resummation near the kinematic threshold limit {xt — 2pr/y^ — > 1), has a 
far smaller sensitivity to scale, compared to NLO calculations, and provides cross sections 
close to those of NLO calculations with a scale of = pr/2. Also, for our energies, the 
calculation of Ref. [67] using /c^-resummation, and of Ref. [63] , which treats simultaneously 
threshold and recoil effects in direct-photon production, yield a substantially larger cross 
section than the NLO result. However, no such calculations are available for inclusive me- 
son production. In their absence, we use a PQCD-based model that incorporates transverse 
kinematics of initial-state partons to study the principal consequences of additional kx for 
high-pT production processes. 

Because the unmodified PQCD cross sections fall rapidly with increasing px, the net 
effect of the "kr smearing" is to increase the expected yield at higher pr's. An exact treat- 
ment of the modified parton kinematics has been implemented in a Monte Carlo calculation 
of the leading-order (LO) cross sections for high-p^ particle production [68], with the kr 
distribution for each of the incoming partons represented by a Gaussian with one adjustable 
parameter ((/cj^)). Unfortunately, no such program is available for NLO calculations, and 
so we approximate the effect of kr smearing by multiplying the NLO cross sections by the 
corresponding LO /cT-enhancement factors. Admittedly, this procedure involves a risk of 
double-counting since some of the /c^-enhancement may already be contained in the NLO 
calculation. However, we expect such double-counting effects to be small. 

The (kr) values used in the calculation of the LO /cT-enhancement factors are similar 
to those employed in comparisons of kinematic distributions in data involving production 
of high-mass 77, 77r°, and 7r°7r° systems with calculations relying on the same LO program 
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(see Refs. [4,5] for further details). For these comparisons, we used the LO versions of 
the CTEQ4 distribution and (where appropriate) BKK fragmentation functions, and an 
average transverse momentum of 0.6 GeV/c for the 7r° mesons relative to the fragmenting 
parton direction (varying this parameter in the range 0.3 to 0.7 GeV/c does not affect our 
conclusions) [69,70]. 

Comparisons of the /cr-enhanced calculations with data at 530 GeV/c are displayed in 
Fig. 17, indicating reasonable agreement for the (kr) values chosen. Similar conclusions can 
be drawn from comparisons between calculations and data at 800 GeV/c, as illustrated in 
Fig. 18. It is interesting to compare the fractional differences between data and the kr- 
enhanced NLO calculations using BKK and KKP fragmentation functions on a linear scale. 
Such comparisons are shown in Fig. 19 for pBe interactions at 530 and 800 GeV/c. The 
/cT-enhanced calculations using the KKP fragmentation functions are seen to reproduce the 
shape of the 7r° cross section better than calculations using the earlier BKK fragmentation 
functions, but small discrepancies are still present. Similar conclusions can be drawn from 
comparisons of NLO QCD with our data on pp collisions (Fig. 20), although in this case 
the fractional differences using the KKP fragmentation functions are systematically greater 
than zero when we employ the same (kr) values used in Fig. 19. It is, however, worth noting 
that the hydrogen nucleus is not typical, having an up to down quark ratio of two rather 
than a value of less than one as in the case of standard nuclei, such as Be (Cu) for which 
u/d — 0.93 (0.94). The theoretical calculations address this difference by calculating the 
nuclear cross sections using the actual numbers of protons and neutrons present in each case, 
but this alone does not fully compensate for the influence of the nuclear media surrounding 
interactions in targets other than hydrogen. As a result, values of (kr) that yield agreement 
with nuclear data adjusted by may not be fully appropriate for hydrogen. In particular, 
raising (kx) by only 0.05 GeV/c drops the points in the lower two quadrants of Fig. 20 by 
enough (fa 0.075 at 4.5 GeV/c) to approximately center these distributions on zero. 

Figures 21 and 22 show the cross sections for inclusive 7r° production as functions of 
rapidity for pBe interactions at 530 and 800 GeV/c, for several intervals in pt- The expected 
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peaking at a scattering angle near 90° in the center of mass (ycm = 0) develops slowly as a 
function of pt- The shapes and normalizations of the data are in good agreement with the 
/cT-enhanced calculations. 

Both theoretical and experimental uncertainties are reduced in the ratio of invariant cross 
sections for 7r° production at 800 and 530 GeV/c, allowing, in principle, a more sensitive 
test of the calculations. Figure 23 displays this ratio compared to conventional ((A;t)=0) 
and /cT-enhanced NLO results using KKP fragmentation functions. The (kx) values used 
here correspond to those used in Figs. 17 and 18. The energy dependence of the data is 
accommodated better by the fcr-enhanced theory. Similar conclusions are obtained for the 
hydrogen target data (not shown). 

The results discussed in this section are not very sensitive to the reasonably well-known 
parton distribution functions [4] (quark distributions being of primary importance here). 
Methods similar to the ones described in this paper have been applied to analyze high-p^ 
hadron spectra inpp, pA, and AA collisions [71,72], and kx effects have been found important 
for describing data on inclusive production of charged mesons. 

B. r) production 

Cross sections for inclusive rj production are tabulated in Tables VIII through XIII. 
Theoretical descriptions of ?7-meson production differ from the tt" case primarily because of 
differences in the fragmentation of partons into the particles of interest. To investigate this 
aspect, we present r^/vr^ relative production rates as functions of pt and ^cm (for two pt 
ranges) in Figs. 24 and 25 — the data average to a value of 0.45 ± 0.01 at 530 GeV/c, and 
0.42 ± 0.01 at 800 GeV/c (for 3 < pr < 8 GeV/c). 

V. SUMMARY 

The invariant cross sections for and 1] production have been measured for pp and 
pBe collisions at 530 and 800 GeV/c as functions of and ycm, over the kinematic range 
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1 < pt < 10 GeV/c, and 1.5 units in rapidity. Results from /cr-enhanced NLO QCD cal- 
culations are in reasonable agreement with our measured 7r° cross sections. Employing the 
recent KKP fragmentation functions in the NLO QCD calculations was found to improve 
the description of the detailed shape of our 7r° cross sections (as a function of pr) relative 
to theoretical results obtained using the earher BKK fragmentation functions. The mea- 
sured rj/i:^ production ratios, which provide information about the relative fragmentation 
of partons into these mesons, are 0.45 ± 0.01 at 530 GeV/c and 0.42 ± 0.01 at 800 GeV/c, 
averaged over pr and rapidity. 
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FIG. 1. Plan view of the Fermilab Meson West spectrometer, as configured for the 1991-92 
fixed target run. 
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FIG. 2. Distribution of the reconstructed vertex position (along Z) for the combined 530 and 
800 GeV/c data. The clearly-resolved individual target elements are, from the left: two copper 
foils, the liquid hydrogen target (a mylar flask enclosed between two thin Be windows), followed 
by the main beryllium target. The two silicon strip detectors (SSDs) immediately upstream and 
four downstream of the target can also be easily distinguished. 
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FIG. 3. The physical layout of the liquid argon electromagnetic calorimeter (EMLAC). 
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FIG. 4. 77 mass distributions in the region of the tt^ (top row) and the r] (bottom row) mesons 
from pBe interactions at 800 GeV/c, for several ranges of 77 pT- values. 
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FIG. 5. pt spectra of tt^ candidates from 530 GeV/c proton interactions on Be, selected via 
several different triggers implemented with significantly different prescale factors. The data have 
only been corrected for trigger prescale factors. 
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FIG. 6. Effect of beam halo muon induced shower discriminators on the 77 invariant mass 
distribution in the vr'^-mass region, for 77 pairs with pT > 7 GeV/c in the 530 GeV/c data (left) 
and in the 800 GeV/c data (right). 
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FIG. 7. Comparisons of distributions in photon Er — E^, for eight photon-energy bins, between 
data (histogram) and the Monte Carlo simulation (points) for the 530 GeV/c sample. The photons 
are from 7r° candidates with pT > 3.5 GeV/c. 
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FIG. 8. Comparison of the photon £^front/^^total distributions, in nine photon-energy bins, 
between data (histogram) and the Monte Carlo simulation (points), for the 800 GeV/c sample. 
The photons originated from decays of 7r° candidates with px > 3.5 GeV/c. 
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FIG. 9. 77 mass distributions in the 7r° signal region in the 530 GeV/c data (histogram) 
compared to Monte Carlo simulations (points) for two requirements on minimum px- 
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FIG. 10. 77 mass distributions in the rj signal region in the 800 GeV/c data (histogram) 
compared to Monte Carlo simulations (points) for two requirements on minimum px- 
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FIG. 11. Comparison of 77 energy asymmetry distribution for tt^ mesons in the data (his- 
togram) and the Monte Carlo (points) for 800 GeV/c pBe interactions, for the pr intervals 
4.0 < pt < 5.5 GeV/c and 5.5 < pT < 7.0 GeV/c. These distributions have been corrected 
for contributions from background sources. 
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FIG. 12. Probability for accepting a candidate (in the 530 GeV/c Monte Carlo sample) as a 
function of pr for three rapidity intervals and different triggers. This probability reflects the losses 
due to the reconstruction algorithm, the -Efront /-E^total criterion, and the individual trigger 
requirements (see text for details). 
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FIG. 13. Invariant differential cross sections (per nucleon) for tt" production as a function of 
7r° pT in pp and pBe interactions at 530 GeV/c. Cross sections have been averaged over the full 
rapidity range, —0.75 < ycm < 0.75. The error bars have experimental statistical and systematic 
uncertainties added in quadrature. 
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FIG. 14. Invariant differential cross sections (per nucleon) for tt" production as a function of 
7r° pT in pp and pBe interactions at 800 GeV/c. Cross sections have been averaged over the full 
rapidity range, —1.0 < ycm < 0.5. The error bars have experimental statistical and systematic 
uncertainties added in quadrature. 



42 



10* 



10 



10 



10 



production at 800 GeV/c 

-i-o<ycm<o-5 

• pBe [pb/(GeV/c) per nucleon] 

□ pp [nb/(GeV/c)^ ] 
Stat and sys uncertainties combined 



10 



10 



10 



10 



-4 



10 



10 



-5 



NLO Theory 
\i = Pj/2 

1^ = Pt 

|i = 2p^ 

BKK fragmentation 

. . . . I . I . I I I . 
4 5 



9 10 11 

p^ (GeV/c) 



FIG. 15. Invariant differential cross sections (per nucleon) for vr*^ production as a function of 
7r° pT in pBe and pp interactions at 800 GeV/c compared to NLO PQCD calculations, with scale 
choices of /Lt = pt/2, Pt, and 2pT- The error bars have experimental statistical and systematic 
uncertainties added in quadrature. 
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FIG. 16. Invariant differential cross sections (per nucleon) for tt" production as a function of 
7r° pT in pBe and pp interactions at 530 GeV/c compared to NLO PQCD calculations, with scale /x 
= pt/2 and BKK and KKP fragmentation functions. The error bars have experimental statistical 
and systematic uncertainties added in quadrature. 
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FIG. 17. Invariant differential cross sections (per nucleon) for tt" production as a function of 
7r° pt in pBe and pp interactions at 530 GeV/c compared to A;T-enhanced NLO PQCD calculations 
with scale fi = pt/2. Comparisons are shown for both BKK and KKP fragmentation functions. 
The error bars have experimental statistical and systematic uncertainties added in quadrature. 
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FIG. 18. Invariant differential cross sections (per nucleon) for vr'' production as a function of 
Pt in pBe and pp interactions at 800 GeV/c compared to fcr-enhanced NLO PQCD calculations, 
with scale ^ = pt/2. Comparisons are shown for both BKK and KKP fragmentation functions. 
The error bars have experimental statistical and systematic uncertainties added in quadrature. 
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FIG. 19. Fractional difference between data and /cr-enhanced NLO PQCD calculations for tt 
production in pBe interactions at 530 and 800 GeV/c as a function of 7r° px- The error bars 
represent only statistical contributions. 
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FIG. 20. Fractional difference between data and fc-r-enhanced NLO PQCD results for tt pro- 
duction in pp interactions at 530 and 800 GeV/c as a function of 7r° pr- The error bars represent 
only statistical contributions. 
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FIG. 21. Invariant cross sections per nucleon for tt production in pBe interactions at 
530 GeV/c. Cross sections are shown versus ?/cm for several intervals in pr- The curves repre- 
sent the fcT-enhanced NLO QCD calculations for (A;r)=1.0 GeV/c and ^ = prf^, using KKP 
fragmentation functions. The error bars have experimental statistical and systematic uncertainties 
added in quadrature. 
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FIG. 22. Invariant cross section per nucleon for production inpBe interactions at 800 GeV/c. 
Cross sections are shown versus ycm for several intervals in pr- The curves represent the 
/cT-enhanced NLO QCD calculations for (A;t)=1.1 GeV/c and /j, = pr/2, using KKP fragmen- 
tation functions. The error bars have experimental statistical and systematic uncertainties added 
in quadrature. 
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FIG. 23. Ratio of invariant cross sections for vr production in pBc interactions at 800 and 
530 GeV/c as a function of tt'^ pr, compared to conventional and A;r-enhanced NLO QCD results 
using KKP fragmentation functions. The error bars reflect only the statistical uncertainties. 
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FIG. 24. Ratios of t] to vr'' invariant cross sections from 530 GeV/c]3Be interactions, as functions 
of (top), and of t/cm (bottom) for two pT ranges. The error bars reflect only the statistical 
uncertainties. 
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FIG. 25. Ratios of t] to vr'' invariant cross sections from 800 GeV/c]3Be interactions, as functions 
of (top), and of t/cm (bottom) for two pT ranges. The error bars reflect only the statistical 
uncertainties. 
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TABLES 



TABLE I. Compilation of systematic uncertainties for inclusive 7r° and r) cross sections. The 
ranges correspond to variations between low-pT to high-px values. 









CiT T Tti r'P>T'1' 511 ■n'^^7" 








o /o 


O /u 




5% 9% 


5% 9% 


Beam Contamination (530 GeV/c beam) 


0% ^ 7% 


0% ^ 7% 


Reconstruction Efficiency 


5% 


6% 


Trigger 


2% ^ 0% 


5% ^ 0% 


Beam Halo Muon Rejection 


1.5% 


1% 


Geometrical Acceptance 


1% 


1% 


Photon Conversions 


1% 


1% 


Vertex Finding 


1% 


1% 


Background Subtraction 


2% ^ 0% 


2% 


Total 


11% ^ 15% 


13% ^ 15% 
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TABLE II. Invariant differential 


cross sections (EcPa/dp^ 


) per nuclcon for the inclusive reac- 


tion pBe 


7r°X at 530 and 800 GeV/c, averaged over the ra 


,pidity intervals —0.75 < Vcrn < 0.75 


and -1.0 < 


Vcm. < 0.5, respectively. 






Pt 




pjje ^ OOU \jr€,V J C 


j^Be @ 800 GeV/c 


(GeV/c) 




[po/[Lfev /c) ) 




1.00 - 1.20 






4- zlQ -k fil 

o'-i.o Zxz ly ^ ux 


1.20 - 1.40 




inn _L 14 _f 21 


^p.'j 4- 90 _L on 
zu i zc zo in ou 


1.40 - 1.60 




(^Q4-|_7Q-|_79 


1 07 _|_ 1 1 -U 1 9 

-LU / IL -L J. in ±z 


1.60 - 1.80 




zu. ( ZJZ 0.0 m Zr.o 


49.3 ± 5.7 ± 5.5 


1.80 - 2.00 




10. -L ZjZ Z._L ZlZ -L-U 


21.5 ± 3.2 ± 2.4 


2.00 - 2.20 




«7_I_-| Q_i_n7i 

•ij. 1 zc L .0 z^ yj. 1 1- 


7.3 ± 1.7 ± 0.82 


2.20 - 2.30 




9 qn _|_ n M _L n 01 


3.92 ± 0.24 ± 0.44 


2.30 - 2.40 






2.67 ± 0.22 ± 0.30 


2.40 - 2.50 






1.77 ± 0.13 ± 0.20 


2.50 - 2.60 




n QAfi -\- n ns^ -1- n 1 n 


1.208 ± 0.094 ± 0.13 


2.60 - 2.70 




U.OOy !E U.U-LO Zl U.UOo 


0.818 ± 0.046 ± 0.091 








inhliGeV lc\^\ 


2.70 - 2.80 




4iD. ( iL 0.0 zL 44 


K4Q _|_ OQ _|_ 79 

uiy 00 ZL 1 ^ 


2.80 - 2.90 




ZoD.o zL 0.0 iL Ol 


zn lu Zl ou 


2.90 - 3.00 




10c; C _L /I _L 01 


ox 1 .U 1 .0 ZL OU 


3.00 - 3.10 




IQK^iJ-Q 1 1 1 4 
ioO.O it o.i it 14 


991 Q _|_ C Q _|_ 9c; 

zzx.o ^ 0.0 ^ ZU 


3.10 - 3.20 




nv A _i_ T 1 in 

y ( .4 it z. i it lu 


ifil 7_|_4c;_LiQ 


3.20 - 3.30 




/ z.y it Z.O it /.O 


11f;9_|_QQ_L1'3 


3.30 - 3.40 




47. y ± 1.4 ± 0.1 


ni f;_|_Q9_Lin 
yi.u ^ o.z ^ ±u 


3.40 - 3.50 




Qfi Q _L 1 rr 1 Q 
00. it 1.0 it 0.0 


f;9 1 _|_ 9 9 _L Q 

UZ._L IC Z.Z ^ U.O 


3.50 - 3.60 




Of; Qn _L n nc _l '7 

zo.oy zt u.yo it z. / 


40 n _|_ 1 _L 4 Q 
10. y ZIZ _L,o 1.0 


3.60 - 3.70 




1 7 QQ _|_ n 01 _L 1 Q 




3.70 - 3.80 




1 Q Qi _|_ n cn + 1 4 


9ft 4 _|_ 1 9 _L 9 q 
zu.t: IL ±.z ID z.y 


3.80 - 3.90 




1 n 1 7 _|_ n 9q J_ 1 1 
±u. J- 1 zn u. zy HZ j. . _l 


90 -1- 1 n -k 9 

ZU.U ^ J-.U ^ Z.O 


3.90 - 4.00 




7 Q4q 4. n r\r:r: n 7K 
1 .oiy Z3Z U-Utjij ziz u. 1 


14.56 lb 0.76 lb 1.6 


4.00 - 4.10 




zlfT^ -1- n flAI -k n 

O.^UO ^ U.U'il ^ U.UO 


10.86 lb 0.44 lb 1.2 


4.10 - 4.20 






8.59 lb 0.12 lb 0.95 


4.20 - 4.30 




O.UOU Z^ U.UZU ZL u.oo 


6.443 lb 0.096 ib 0.71 


4.30 - 4.40 




9 9^4 -1- n 091 -1- n 9A 

Z.Z04 ^ U.UZX ZL U.Z^ 


4.706 lb 0.073 ib 0.52 


4.40 - 4.50 




1 -i-nni7-i-n is 


3.657 lb 0.067 lb 0.41 


4.50 - 4.60 




1 97Q -un niA-i-n ^A 


2.831 lb 0.054 ib 0.32 


4.60 - 4.70 




n qfiq -un ni9-l-n ID 


2.159 lb 0.046 ib 0.24 






(r\h 1 ( Ht^M 1 r-\'^\ 

[po/yijev /c) ) 




4.70 - 4.80 




'70 A _L in _L '7n 

io4 zt lu ± /y 


I 797 _|_ Q7 _L 1 nn 

I I z 1 ^ o 1 in lyu 


4.80 - 4.90 




DOD.O ± 0.0 ± bU 


1 Q4Q _|_ 94 _L 1 c;n 


4.90 - 5.00 




/lOT" Q _L "7 1 _L /1'7 
4z;.oit (.llt4/ 


1 n4q H- 99 + 1 90 


5.00 - 5.10 




000. U ± D.O ± o7 


QI 1 _|_ 1 C _L qi 

oil ^ lo ^ yi 


5.10 - 5.20 




oc^/i n _L _L OQ 
Z04.U it O.z it zo 


ftcro _|_ 1 7 _L 7/1 


5.20 - 5.30 




1 n7 Q _L /I Q _L 00 

ly /.o it 4.0 it ZJ. 


4Q« _|_ 1 « _L ccp; 

lOU IL lU HI OO 


5.30 - 5.40 




1 Kn K _L Q n 1 1 '7 
loU.O it v5.y it 1 / 


AAA -|_ 1 Q 4- A(\ 

4: 1 1 31 1 HI "410 


5.40 - 5.50 




ii7n-i-'i/i-i-i'i 
11/. U it v5.4 it lv5 


Ol 9 _L 1 n _L or: 
o 1 z HZ 1 u zn o o 


5.50 - 5.60 




00. D it z.y zt y.o 


9/1 q 7 _L 7 Q 4- 90 

Z-^IO. I ZL i .0 ZtZ Zo 


5.60 - 5.70 




fiO /I J_ rr 1 n 

DZ.4 zt z.o zt D.y 


907 1 _|_ ft Q _L 94 

ZU 1 .1 ^ U.O ^ Z4: 


5.70 - 5.80 




00.0 ± Z.O ± O.Z 


lc;qr:_|_fti_LiQ 
A-Oa.O ^ U.l ^ 10 


5.80 - 5.90 




43.5 lb 2.0 lb 4.9 


l9Sfi+'^fl+lc; 

-Li;0.1J _L O.VJ 1 XiJ 


5.90 - 6.00 




35.1 lb 1.8 lb 3.9 


-HJO.^ _L 0.0 ' Xii 


6.00 - 6.25 




22.42 lb 0.85 lb 2.5 


74 ^4-904-8 ^ 
Cri.o z.u ZL o.tj 


6.25 - 6.50 




11.98 lb 0.61 lb 1.4 


43.9 ±1.6 ±5.1 


6.50 - 6.75 




5.84 lb 0.41 lb 0.67 


23.5 ±1.1 ±2.7 


6.75 - 7.00 




4.32 lb 0.35 lb 0.50 


15.58 ± 0.83 ± 1.8 


7.00 - 7.50 




1.65 lb 0.15 lb 0.20 


7.31 ± 0.40 ± 0.87 


7.50 - 8.00 




0.336 lb 0.066 lb 0.041 


2.56 ±0.25 ±0.31 


8.00 - 9.00 




0.099 zt 0.024 ±0.013 


0.72 ±0.11 ±0.09 


9.00 - 10.00 




0.0077 lb 0.0055 lb 0.0010 


0.068 ± 0.024 ± 0.009 


10.00 - 12.00 






0.020 ± 0.019 ± 0.003 
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TABLE III. Invariant differential cross sections [Ed?a / dp^) for the inclusive reaction pp — > 
7r°X at 530 and 800 GcV/c, averaged over the rapidity intervals —0.75 < ^cm < 0.75 and 



-1.0 < J/cm ^ 0.5, respectively. 



Pt 


pp @ 530 GeV/c 


pp @ 800 GeV/c 


{GeV/c) 


(/.b/{Gey/c)2) 


(M6/(GeV/c)2) 


1.00 - 1.40 


284 ± 44 ± 30 


531 ± 72 ± 59 


1.40 - 1.80 


63 ± 10 ±6.7 


113 ± 17 ± 13 


1.80 - 2.20 


13.5 ±3.6 ±1.4 


17.3 ±5.2 ±1.9 


2.20 - 2.40 


3.36 ± 0.81 ± 0.36 


3.25 ±0.43 ±0.36 


2.40 - 2.60 


1.15 ±0.14 ±0.12 


1..33± 0.22 ±0.15 




{nb/(GeV/c)'-') 


{nb/{GeV/cY) 


2.60 - 2.80 


475 ±17 ±51 


693 ± 57 ± 77 


2.80 - 3.00 


205.2 ± 7.4 ± 22 


315 ± 13 ± 35 


3.00 - 3.20 


98.9 ±4.6 ±11 


171.7 ±7.7 ±19 


3.20 - 3.40 


52.5 ±3.3 ±5.6 


91.7 ±4.7 ±10 


3.40 - 3.60 


25.6 ±1.7 ±2.7 


47.4 ±3.5 ±5.2 


3.60 - 3.80 


12.63 ± 0.99 ± 1.3 


25.7 ±2.3 ±2.8 


3.80 - 4.00 


7.22 ± 0.24 ± 0.77 


14.7 ±1.3 ±1.6 


4.00 - 4.20 


4.170 ±0.060 ±0.45 


9.18 ±0.65 ±1.0 


4.20 - 4.40 


2.263 ± 0.037 ± 0.24 


5.03 ±0.15 ±0.56 


4.40 - 4.60 


1.279 ±0.027 ±0.14 


2.906 ± 0.099 ± 0.32 




(ph/{GeV/cY) 


{ph/{GeV/cY) 


4.60 - 4.80 


757 ± 18 ± 82 


1688 ± 64 ± 190 


4.80 - 5.00 


450 ± 13 ± 49 


1057 ± 36 ± 120 


5.00 - 5.20 


279.6 ±9.5 ±31 


600 ± 27 ± 67 


5.20 - 5.40 


164.3 ±7.6 ±18 


400 ± 21 ± 45 


5.40 - 5.60 


101.4 ±5.4 ±11 


256 ± 17 ± 29 


5.60 - 5.80 


59.9 ±4.0 ±6.7 


163 ± 11 ± 19 


5.80 - 6.00 


38.3 ±3.3 ±4.3 


98.4 ±7.9 ±11 


6.00 - 6.25 


22.5 ±2.1 ±2.5 


75.8 ±5.3 ±8.7 


6.25 - 6.50 


11.1 ±1.5 ±1.3 


34.2 ±3.1 ±4.0 


6.50 - 6.75 


5.53 ± 0.94 ± 0.64 


27.6 ±2.6 ±3.2 


6.75 - 7.00 


4.15 ±0.80 ±0.48 


12.3 ±1.9 ±1.5 


7.00 - 7.50 


1.35 ±0.31 ±0.16 


8.2 ±1.2 ±0.98 


7.50 - 8.00 


0.70 ± 0.23 ± 0.09 


3.35 ± 0.65 ± 0.41 


8.00 - 9.00 


0.185 ±0.093 ±0.023 


0.67 ±0.19 ±0.08 


9.00 - 10.00 




0.083 ± 0.098 ± 0.011 


10.00 - 12.00 




0.016 ± 0.016 ± 0.002 
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TABLE IV. The averaged invariant differential cross section per nucleon {Ed^a/dp^) as a 
function of rapidity and pt for the inclusive reaction pBe tt'^X at 530 GeV/c. 



S/cm 






Pt (LieV/cJ 








1.00 - 1.50 


1.50 - 2.00 


2.00 - 2.50 


2.50 - 3.00 






(/ib/(GeV/cj ) 


(pb/(GeV /c) ) 


(fib/(GeV/c) ) 


[710/ \ij&V / C) ) 


-0.750 - 
-0.625 - 


-0.625 
-0.500 


388 ± 43 ± 41 


27.2 ± 7.9 ±2.9 


6.5± 1.6±0.7 


820 ± 180 ± 90 

AC\Q -\- -\- AA 

4uy it 00 it 44 


-0.500 - 
-0.375 - 


-0.375 
-0.250 


313 ± 40 ±33 


25.5 ± 6.0 ±2.7 


2.9± 1.2±0.3 


490 ± ih 00 
^ xu zl 00 


-0.250 - 
-0.125 - 


-0.125 
0.000 


324 ± 33 ±34 


29. 7± 6.2 ±3.1 


4.1 ± 1.2 ±0.4 


CI q J_ 1 9 _|_ 
010 m xz It 00 

P;87 4-1^-1- 
Oo 1 IC Xo It Do 


0.000 - 

0.125 - 


0.125 

0.250 


233 ± 30 ±25 


33.4 ± 5.7 ±3.5 


5.1 ± 1.2 ±0.5 


CO? Q _!_ Q Q -U c^^; 
oz/.oit y.o it 00 

OUz.U it o.Z it 04 


0.250 - 
0.375 - 


0.375 
0.500 


241 ± 29 ±25 


31.5 ± 4.9 ±3.3 


2.7± 1.3 ±0.3 


^00.0 It (.Dm 4y 

/I OV _!_ 7 PC _|_ /I fl 
4Z(.Olt (.Olt4D 


0.500 - 
0.625 - 


0.625 
0.750 


293 ± 29 ±31 


17.5 ± 5.0 ±1.9 


2.28 ± 0.97 ±0.24 


001. Ult O.U it 41 

Q/11 Q _!_ T _|_ r>(; 
o4i.oit (.oitoD 






o.UU — o.oU 


Q fin A nn 
o.oU — 4.UU 


4.UU — 4.0U 


4.0U — O.UU 






(no/(G-ev/cJ ) 


[nb/yLjeV/c) ) 


[no/ [LreV /c) ) 


(no/ (^CrCV /CJ ) 


-0.750 - 


-0.625 


CO 1 AC 1 '7 1 

Do.o ± 4.0 ± 1.6 


9.00 ± O.DD ± i.O 


i~i COT 1 r\crz 1 r\ OT 

z.ozT ± U.Ubo ± 0.27 


U.ooo± O.UzU ± U.Odo 


-0.625 - 


-0.500 


oo.o III 0.0 III y. i 


1*^0-1- 1 9-1-1 ^ 

10. y m i.z m i.o 


9 oaa -I- n nR9 -1- n "^i 


u.D( / zn. u.uiy it u.uio 


-0.500 - 


-0.375 


7o.o ± 5.1 ± 0.4 


14. z ± l.U ± 1.0 


c\c A 1 n r\co 1 f\ OCT 

3.2o4 ± U.Uoo lb U.oo 


n TOO 1 n ni o i n nocr 


-0.375 - 


-0.250 


Qi n_L. 094.07 


X 1 . X ^ X . X ZLI X .0 


o.oyy ^ u.uox ziz u.oy 


n qnc; -i- n m 4- n riQS 

u.yuo 31 u.uxo Zl u.uyo 


-0.250 - 


-0.125 


89.5 ± 3.3 ±9.6 


17.8 ± 1.1 ±1.9 


3.864 ± 0.042 ±0.41 


0.963 ± 0.017 ±0.10 


-0.125 - 


0.000 


Q1 /1_L.QQ_l_OQ 

yi.4i:it 0.0 it y.o 


on n _i_ 1 _i_ 1 
zU.U it 


A 1 Qc _i_ n H/i 1 -J- n /I 

4.ic50 ± U.U4i it U.40 


1 fiQa 4- n m ? _i- n 1 1 
i.UocS ± u.Uii iru.ii 


0.000 - 


0.125 


88.5 ± 2.9 ±9.4 


17.30 ± 0.96 ±1.8 


4.001 ± 0.038 ± 0.43 


0.927 ± 0.015 ±0.10 


0.125 - 


0.250 


( 1 . 1 ± Z.^ it cS.o 


IK Al -X- n 'TQ _!_ 1 K 
10.4 (itU. (Oltl.D 


0. 1 ZU it U.UoO 31 U.4U 


n QQ1 -1- nnic;_i_n in 
U.yoi it U.UiO it U.iU 


0.250 - 


0.375 


78.9 ± 2.5 ±8.4 


16.77 ± 0.85 ±1.8 


3.420 ± 0.032 ± 0.37 


0.872 ± 0.015 ± 0.094 


0.375 - 


0.500 


( 4.0 it Z.D it o.U 


1 Q 1 C _!_ n '7'7 1 1 A 

lo.lOit U. / ( it i.4 


o.loU it U.Uo4 it U.o4 


n '7A'i -L- n ni /i _i_ n non 
U. ( 4o it U.U14 it U.Uc5U 


0.500 - 


0.625 


65. 1± 2.6 lb 7.0 


13.9 lb 1.0 ±1.5 


2.670 ± 0.033 ±0.29 


0.629 ± 0.014 ±0.068 


0.625 - 


0.750 


00.0 =t z.D zt O.y 


in 1 /I _i_ n Qn 1 1 1 
iU.14it U.oy it 1.1 


1 cn -L n nun -U n oq 
z.iOy it U.UoU ± U.zo 


u.4Dy =t u.uio =t u.uoi 






5.00 - 5.50 


5.50 - 6.50 


6.50 - 8.00 


8.00 - 10.00 






(pb/{GeV /cy) 


(pb/(GeV /c)^) 


[pb/{GeV/c) ) 


(pb/{GeV/c} ) 


— U. ( ou 
-0.625 - 


— U.DZO 

-0.500 


144.3 ± 7.6 ± 16 
202.1 ± 8.6 ± 22 


24.9 ± 2.0 ± 2.8 
28.9 ± 1.9 ±3.2 


1.51 ± 0.22 ±0.18 


0.015 ± 0.027 ±0.002 


-0.500 - 
-0.375 - 


-0.375 
-0.250 


215.9 ± 7.7 ± 24 
236.4 ± 7.8 ± 26 


40.0 ± 2.1 ±4.5 
44.0 ± 2.2 ±4.9 


2.02 ± 0.23 ±0.24 


0.053 ± 0.027 ±0.007 


-0.250 - 
-0.125 - 


-0.125 
0.000 


253.7 ± 7.9 ± 28 
268.0 ± 8.2 ±29 


48.4 ± 2.3 ± 5.4 
50.2 ± 2.3 ±5.6 


3.31 ± 0.30 ±0.39 


0.059 ± 0.030 ± 0.007 


0.000 - 
0.125 - 


0.125 
0.250 


251.5 ± 7.4 ±28 
268.7 ± 7.7 ±29 


46.1 ± 2.1 ±5.2 
45.5 ± 2.1 ±5.1 


3.04 ± 0.27 ±0.36 


0.152 ± 0.050 ± 0.019 


0.250 - 
0.375 - 


0.375 
0.500 


230.7 ± 7.2 ±25 

182.8 ± 6.6 ±20 


42.7 ± 2.0 ±4.8 
34.7 ± 1.9 ±3.9 


2.68 ± 0.28 ±0.32 


0.029 ± 0.021 ± 0.004 


0.500 - 
0.625 - 


0.625 
0.750 


164.0 ± 6.6 ± 18 
111.0 ± 5.7 ± 12 


24.8 ± 1.8 ±2.8 
15.6 ± 1.4 ±1.7 


1.58 ± 0.23 ±0.19 


0.014 ± 0.014 ±0.002 
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TABLE V. The averaged invariant differential cross section [Ed^a /dp^) as a function of rapid- 
ity and pt for the inclusive reaction pp — > tt^X at 530 GeV/c. 



S/cm 






Pt (LieV/cj 










1.00 - 2.50 


2.0U — o.UU 


Q nn Q c^n 
O.UU — o.oU 


Q cn A nn 
O.oU — 4.UU 








yrio / [kjCv f c) ) 






-0.750 - 
-0.625 - 


-0.625 
-0.500 


147 ± 


15 ± 16 


izAn 1 ion 1 cr\ 

54U ± loU ± oU 

Qc^n -J- QQfi -U 1 nn 
you it ooU it lUU 


£;o _J_ 1 n _i_ T 

DO ± lU ih 7 

OD it lo it 1 


0.7 ih i.4 ± U.y 
inn-i-i yi-ui 1 

lU.U it 1.4 it 1.1 


-0.500 - 


-0.375 


119 ± 


15 ± 13 


Qtjc _L /in _L Qn 

ODD ± 4u ± oy 


^1 _L 1 n _i_ 
oi ± iU ± D 


1 O Q -L O '7 -L 1 Q 

iz.o ± Z.7 ± i.o 


-0.375 - 


-0.250 


AAn _j_ QQ _|_ 47 


«c c_|_ 7o_|_7n 


14Q4- 9^-1-1 fk 


-0.250 - 


-0.125 


123 ± 


12 ± 13 


"4:00 It ZD It 


lO. 1 It o.y It C5.Z 


10.41: It Z.'i It l.D 


-0.125 - 


0.000 


Oo^ It di 1 


Ql 7_|_70_|_Q7 
C51. 1 It < -Z It o. ( 


i^.o It z.o di i.o 


0.000 - 


0.125 


91 ± 


11 ±10 


A^T -L- on -U AQ 
40 / it Z\J it 4y 


1 o. 1 it O.O it 0.4 


IQCI-UOO-UI A 
lo.Olt Z.Z It 1.4 


0.125 - 


0.250 


A 0/1 -1- 1 O -1- /I 


K7 Q 4- /I -1- T' O 
D( .o it 0.4 zt i .Z 


lo.O it l. ( it 1.4 


0.250 - 


0.375 


96 ± 


11 ±10 


QQQ J_ 1 Q -U /19 

oyo It io It ^tz 


Do.O It 0. ( It D.O 


11 9-1- 1 j^-l-l 9 


0.375 - 


0.500 


Q(^0 _|_ 1 7 _|_ QO 

vSoy it i ( it oy 


PS) 7-1- c: O -L 7 
Dz. (it O.y it D. / 


i/ii-i-9n-i-i 
14.1 it z.U it 1.0 


0.500 - 


0.625 


116 ± 


11 ± 12 


04/ it lo it O / 


f^7 O -U 1 -U 1 
O/.Z it 0.1 it 0.1 


11 1J-99-I-1 O 
11.1 It Z.Z It l.Z 


0.625 - 


0.750 


QnQ _l_ 1 Q _l_ QO 
oUo it lo it OZ 


AO, A A- O -U Q 

4y.4 it O.y it O.O 


( .o it 1.0 it U.o 






4.00 - 4.50 


A c;n nn 
4.0U — O.UU 


O.UU — O.OU 


c;n fi f^n 
O.oU — O.OU 






(nb/{GeV/cf) 


[no/ [Lrev / c) ) 


{pb/{GeVlc] ) 


(p6/(GeV/c) ) 


-0.750 - 


-0.625 


2.02± 


0.15 ±0.22 


U.blU ± U.Udz ih U.UoD 


iU4 ± io ziz li 


z7.o ± 4.0 ± 6.1 


-0.625 - 


-0.500 


2.46 ± 


0.14 ±0.26 


U-010 m U.U4tD It U.UOD 


1 _|_ 9 1 J- 1 Q 
lOD zt Zl It lo 


Z±.U ^ '4:.Z ^ Z.4t 


-0.500 - 


-0.375 


2.74 ± 


0.11 ±0.29 


U.o45 lb U.UoD ± U.U7U 


on'? 1 1 o 1 oo 

207 ± 18 ± 23 


40.5 ± 5.4 ± 4.6 


-0.375 - 


-0.250 


3.34 ± 


0.12 ±0.36 


n QQQ _j_ r\ n^c 4- n nQfl 


1 no -1- 1 7 -1- OI 




-0.250 - 


-0.125 


3.318 ± 


0.091 ± 0.36 


0.809 ± 0.037 ± 0.087 


231 ± 18 ± 25 


41.0 ± 5.2 ±4.6 


-0.125 - 


0.000 


3.555 ± 


0.095 ± 0.38 


n oi ? _L n n/in _l n nno 

u.yii it U.U4U it u.uyy 


OQO _i- 1 n -1- Of? 

zoy ± ly ± zo 


^io n 1 -1 1 T 1 
Oz.y it D.4 it / .1 


0.000 - 


0.125 


3.565 ± 


0.092 ± 0.38 


0.807 lb 0.034 lb 0.087 


266 lb i9 lb 29 


52.iib 5.5 lb 5.8 


0.125 - 


0.250 


3.316 ± 


0.083 ± 0.36 


n _L n nQT _l n nno 
U.o4o ± U.UoY ± U.Uyz 


O/in -L 1 Q _J- 0'7 

z4y ± io ± 27 


CO c -L c a -L_ c n 
oz.o ± o.D ± o.y 


0.250 - 


0.375 


2.960 ± 


0.072 ± 0.32 


0.723 lb 0.033 lb 0.078 


235 lb 18 ± 26 


38.7 lb 4.8 lb 4.3 


0.375 - 


0.500 


2.791 ± 


0.078 ± 0.30 


n fiQQ _L n nQQ _l n i\'7A 
U.ooo it U.Uoo it U.U ( 4 


onn J- 1 '7 _i_ oo 
zUU it 1 ( it zz 


on T 1 'I O 1 Q Q 

zy. / it 4.Z it O.O 


0.500 - 


0.625 


2.356 ± 


0.075 ± 0.25 


0.624 ± 0.034 lb 0.067 


i79 zt i7 zb 20 


25.3 lb 4.2 ±2.8 


0.625 - 


0.750 


1.917± 


0.069 ± 0.21 


n /ic;i _L n nQn _L n n/in 
U.4ol ± U.UoU ± U.U4y 


11/^ _L 1 A -L "1 Q 

iio ± i4 ± io 


OI Q -L A 1 -L O /I 
zi.o it 4.1 it Z.4 






6.50 - 8.00 












(p6/(GeV/c)2) 








— U. ( oU — 
-0.625 - 


— U.OzO 

-0.500 


2.00 ± 


0.57 ± 0.24 








-0.500 - 
-0.375 - 


-0.375 
-0.250 


2.62 ± 


0.66 ± 0.31 








-0.250 - 
-0.125 - 


-0.125 
0.000 


3.00 ± 


0.70 ± 0.35 








0.000 - 
0.125 - 


0.125 
0.250 


3.62 ± 


0.70 ±0.43 








0.250 - 
0.375 - 


0.375 
0.500 


1.70 ± 


0.50 ±0.20 








0.500 - 
0.625 - 


0.625 
0.750 


0.91 ± 


0.41 ±0.11 
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TABLE VI. The averaged invariant differential cross sections per nucleon {Ed^a/dp^) as a 
function of rapidity and pt for the inclusive reaction pBe vr'^X at 800 GeV/c. 



S/cm 






Pt {KjG\/c) 








1.00 - 1.50 


1.50 - 2.00 


2.00 - 2.50 


2.50 - 3.00 








W^KGeV /c) ) 


(Hb/(GeV/c) ) 


[TlOj [ijCV j C) ) 


-1.00 - 
-0.875 - 


-0.875 
-0.750 


444 ± 71 ± 49 


36 ± 11 ±4 


6.5± 1.9 ±0.7 


477 ± 87 ±53 

QAn -I- 1 an J- on 
o4U it icSU it yu 


-0.750 - 
-0.625 - 


-0.625 
-0.500 


297 ± 65 ±33 


58.8 ± 9.0 ±6.5 


5.9 ± 1.8 ±0.7 


til a _L cc; _L CO 
DlO ± 00 ± Do 

ftoft _|_ cin 4- 7f; 

UOU ^ OU ^ 1 \J 


-0.500 - 
-0.375 - 


-0.375 
-0.250 


372 ± 59 ±41 


57.1 ± 9.4 ±6.3 


1.3± 1.8 ±0.1 


am -I- c;n 4- fic 
DIU it OU it Do 

/ UO 3: 01 IE 1 


-0.250 - 
-0.125 - 


-0.125 

0.000 


304 ± 53 ±34 


55.1 ± 8.5 ±6.1 


4.1 ± 1.6 ±0.5 


4- AK -\- 91A 
1 Oo it 40 it o4 

7QO -1- /Ift -1- Q1 
i It 40 It ol 


0.000 - 
0.125 - 


0.125 
0.250 


269 ± 46 ±30 


45.3 ± 7.7 ±5.0 


5.4 ± 1.6 ±0.6 


1 DO zL ^1 ztz oO 

7Qn 4- 1 nn 4- cn 
* oU it lUU it oU 


0.250 - 
0.375 - 


0.375 
0.500 


409 ± 45 ±45 


31.6 ± 7.6 ±3.5 


4.5 ± 1.4 ±0.5 


i OU it 4y it oo 
f^Qfi -I- /10 -I- ftp; 

Ooo it 4Z it DO 






o.UU — v>.OU 


o.oU — 4.UU 


A nn A c;n 
4.UU — 4.0U 


/I c;n c; nn 
4.0U — O.UU 






{noj {yjiG-v 1 c) ) 


[no/ (CtC V / c) ) 


(nb/[GeV/c) ) 


[no/ [Lrev / c) ) 


-1.00 - 


-0.875 


T r\r\ cr 1 o t i i t 


on n 1 n 1 


/ini 1 noo 1 r\ A A 

4.U1 ± yj.2Q ± U.44 


1 non 1 n n o / ■ i n i o 


-0.875 - 


-0.750 


19zl 1-1- 7 « -\- A A 


99 4 -k 9 9 -k 9 


A 1'\ 4- fl '^7 -1- n f^9 
"4:. / It U.Ol It U.OZ 


"1 1 on _|_ n nQft _|_ n i q 
l.loUlt U.UoOlLU.lO 


-0.750 - 


-0.625 


102.3 ± D.l ± 11 


cr 1 cr 1 00 

2o.6 ± 2.0 ± z.o 


c cr 1 n 01 1 n co 

0.58 ± 0.31 ± O.oi 


1 CIO 1 nnoo i niT 

1.51.2 ± U.Uoo ± U.17 


-0.625 - 


-0.500 


1104-1- co-i-iq 


oco-L 1 04.90 




1 7QC J_ fl n^o _|_ n on 


-0.500 - 


-0.375 


143.0 ± 6.4 ± 16 


30.0 ± 2.1 ±3.3 


7.40 ± 0.40 ± 0.82 


2.029 ± 0.037 ±0.23 


-0.375 - 


-0.250 


1441:. it 0. 1 it ID 


QO Q _1_ 1 _|_ Q C 

oz.o it Z.L it o.D 


o.UU it U.oo it U.oy 


o 1 ni 4- n nQc^ 4- n oq 
z.iui ± U.Uoo it U.zo 


-0.250 - 


-0.125 


149.0 ± 5.7 ± 17 


28.1 ± 2.2 ±3.1 


7.91 ± 0.32 ± 0.88 


2.168 ± 0.036 ±0.24 


-0.125 - 


0.000 


i4l:4t.4l: It O.CS It ID 


zy.y it z.o it 0.0 


i.iOlt U.O( it U.oD 


o 1 Qo 4- n nQQ 4- n o/i 
z.iyz it U.Uoo it U.Z4I: 


0.000 - 


0.125 


142.1 ± 4.8 ±16 


32.1 ± 1.8 ±3.5 


7.96 ± 0.29 ± 0.88 


2.129 ± 0.087 ±0.24 


0.125 - 


0.250 


14t5.4 It -rt-D It iD 


OQ _1_ 1 V _1_ Q Q 

zy.O it 1. ( it 0.0 


T fio J- n OQ J- n Q/1 

/.Dz it U.zo it U.o4 


no 4- n 1 n 4- n oq 
z.\jz it u.iu it yj.zo 


0.250 - 


0.375 


129.3 ± 4.9 ± 14 


29.5 ± 1.8 ±3.2 


6.88 ± 0.22 ±0.76 


1.877 ± 0.086 ± 0.21 


0.375 - 


0.500 


1 OT Q -U IT fr 1 1 A 
LA i .0 ± 0.0 it 14 


on J_ 1 1 Q 
zy.z it z.l it o.Z 


D./Dlt U.Z4ztU. (O 


1 coc _L n nQ/1 -1- n on 
l.oZo it U.UoA ± U.zU 






5.00 - 5.50 


5.50 - 6.50 


6.50 - 8.00 


8.00 - 10.00 






(pb/(GeV/c) ) 


(pb/(GeV/c)^) 


(pb/{GeV/c) ) 


{pb/(GeV/c) ) 


— i.uu — 
-0.875 - 


— U.O 1 o 

-0.750 


268 ± 25 ± 30 
326 ± 16 ±37 


54.3 ± 3.7 ±6.2 
68.1 ± 4.3 ±7.8 


3.54 ± 0.42 ±0.42 




-0.750 - 
-0.625 - 


-0.625 
-0.500 


411 ± 16 ± 46 
505 ± 19 ± 57 


82.8 ± 4.1 ±9.5 
109.3 ± 5.6 ±13 


7.51 ± 0.63 ±0.89 


0.10 ± 0.10 ±0.01 


-0.500 - 
-0.375 - 


-0.375 
-0.250 


595 ± 18 ± 67 
635 ± 17± 72 


120.5 ± 4.9 ±14 
141.5 ± 5.3 ±16 


11.18 ± 0.70 ± 1.3 


0.50 ± 0.11 ±0.06 


-0.250 - 
-0.125 - 


-0.125 
0.000 


646 ± 18 ±73 
639 ± 19 ±72 


142.2 ± 5.0 ±16 
135.2 ± 5.2 ±15 


12.70 ± 0.75 ± 1.5 


0.60 ± 0.16 ±0.07 


0.000 - 
0.125 - 


0.125 
0.250 


645 ± 28 ± 73 

646 ± 37 ±73 


137.7 ± 6.1 ± 16 
132.4 ± 7.2 ±15 


13.96 ± 0.78 ± 1.7 


0.76 ± 0.15 ±0.09 


0.250 - 
0.375 - 


0.375 
0.500 


583 ± 29 ±66 
526 ± 27 ±59 


134.1 ± 6.5 ±15 
107.0 ± 6.0 ±12 


9.93 ± 0.73 ±1.2 


0.48 ± 0.16 ±0.06 
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TABLE VII. The averaged invariant differential cross sections [Ed^a/dp^) as a function of 
rapidity and pt for the inclusive reaction pp — tt^X at 800 GeV/c. 



Vcm PT (GeV/c) 





1.00 - 2.50 


2.50 - 3.00 


3.00 - 3.50 


3.50 - 4.00 




(Mb/(Gey/c)2) 


(nb/{GeV/cf) 


inb/{GeV/cf) 


{nb/{GeV/cf) 


-1.00 - -0.875 


166 ± 


26± 18 


980 ± 580 ±110 


83± 16 ±9 


19.2 ± 7 A ±2.1 


-0.875 - -0.750 


410 ± 110 ±50 


108 ± 14± 12 


18.7± 4.1 ±2.1 


-0.750 - -0.625 


120 ± 


24± 13 


420 ± 63 ± 47 


101 ± 14± 11 


36.9 ± 6.9 ±4.1 


-0.625 - -0.500 


570 ± 100 ± 60 


111± 12± 12 


18.0 ± 3.6 ±2.0 


-0.500 - -0.375 


146 ± 


21 ± 16 


560 ± 110 ± 60 


144 ± 14 ± 16 


26.6 ± 5.5 ±2.9 


-0.375 - -0.250 


600 ± 99 ± 67 


130 ± 13 ±14 


37.1 ± 6.0 ±4.1 


-0.250 - -0.125 


130 it 


20± 14 


760 ± 170 ± 90 


117± 14 ±13 


31.0± 5.2 ±3.4 


-0.125 - 0.000 


640 ± 110 ±70 


117 ± 12± 13 


20.6 ± 4.1 ±2.3 


0.000 - 0.125 


115 ± 


17± 13 


669 ± 95 ± 74 


127 ± 11 ± 14 


23.7 ± 4.5 ±2.6 


0.125 - 0.250 


720 ± 110 ±80 


130 ± 10 ± 14 


26.3 ± 4.0 ±2.9 


0.250 - 0.375 


164 ± 


16 ± 18 


511 ± 84 ±57 


115 ± 11 ± 13 


22.5 ± 3.6 ±2.5 


0.375 - 0.500 


699 ± 97 ± 78 


102 ± 12± 11 


22.4 ± 5.1 ±2.5 




4.00 - 4.50 


4.50 - 5.00 


5.00 - 5.50 


5.50 - 6.50 




(n6/(Gey/c)2) 


{nb/(GeV/c)'^) 


{pb/{GeV/cf) 


(pV{GeV/c)2) 


-1.00 - -0.875 


3.02 ± 


0.25 ± 0.34 


1.18 ± 0.12 ±0.13 


280 ± 46 ± 32 


30.3 ± 5.2 ±3.5 


-0.875 - -0.750 


4.66 ± 


0.65 ± 0.52 


0.957 ± 0.072 ±0.11 


217 ± 25 ± 24 


70 ± 13 ± 8 


-0.750 - -0.625 


8.4 ± 


2.2 ±0.9 


1.350 ± 0.085 ±0.15 


331 ± 34 ± 37 


60.2 ± 8.2 ±6.9 


-0.625 - -0.500 


7.4 ± 


1.5 ±0.8 


1.504 ± 0.093 ±0.17 


423 ± 36 ±48 


99± 12± 11 


-0.500 - -0.375 


6.41 ± 


0.63 ± 0.71 


1.808 ± 0.087 ±0.20 


479 ±41 ±54 


84.9 ± 9.7 ±9.7 


-0.375 - -0.250 


6.86 ± 


0.61 ± 0.76 


1.901 ± 0.081 ±0.21 


517 ± 39 ±58 


124 ± 12 ±14 


-0.250 - -0.125 


6.00 ± 


0.20 ± 0.66 


1.963 ± 0.083 ±0.22 


595 ± 41 ± 67 


110 ± 11 ± 13 


-0.125 - 0.000 


6.14 ± 


0.35 ± 0.68 


1.878 ± 0.081 ±0.21 


552 ± 38 ± 62 


140 ± 13± 16 


0.000 - 0.125 


8.2 ± 


1.0 ±0.9 


1.49 ± 0.13 ±0.17 


607 ± 68 ±68 


150 ± 17 ±17 


0.125 - 0.250 


6.19 ± 


0.57 ±0.69 


1.71 ± 0.22 ±0.19 


567 ± 88 ±64 


111± 14 ±13 


0.250 - 0.375 


7.08 ± 


0.57 ±0.78 


1.68 ± 0.20 ± 0.19 


459 ± 61 ± 52 


97± 12± 11 


0.375 - 0.500 


5.69 ± 


0.52 ± 0.63 


1.81 ± 0.21 ±0.20 


502 ± 64 ± 57 


129 ± 16± 15 




6.50 - 8.00 


8.00 - 10.00 








(pb/iGeV/c)^) 


(p6/(Gey/c)2) 






1 nn n 87^; 

— l.UU — U.OiO 

-0.875 - -0.750 


4.20 ± 1.00 ±0.50 


0.085 ± 0.085 ± 0.010 






-0.750 - -0.625 
-0.625 - -0.500 


8.8 ± 


1.6 ±1.0 


0.23 ± 0.17 ±0.03 






-0.500 - -0.375 
-0.375 - -0.250 


14.0 ± 


1.8 ± 1.7 


0.53 ± 0.24 ±0.07 






-0.250 - -0.125 
-0.125 - 0.000 


13.8 ± 


1.7± 1.6 


0.95 ± 0.43 ±0.12 






0.000 - 0.125 
0.125 - 0.250 


12.7 ± 


2.1 ± 1.5 


0.66 ± 0.30 ±0.08 






0.250 - 0.375 
0.375 - 0.500 


9.6 ± 


1.7± 1.1 
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TABLE VIII. Invariant differential cross sections (Ed^a/dp^) per nucleon 


for the inclusive 


reaction pBe 


r/X at 530 and 800 GeV/c, averaged over the rapidity intervals — 


0.75 < < 0.75 


and —1.0 < Vrm 


< 0.5, respectively. 




Pt 


pBe ® 530 GeV/c 


pBe @ 800 GeV/c 


(GeV/c) 


{nh/(GeV/cf) 


(nh/(GeV/r')'^) 


3.00 - 3.20 


49.2 ±9.0 ±6.0 


QQ _|_ 1 Q _L 1 1 


3.20 - 3.40 


24.4 ±4.5 ±2.9 


'\A 7 A- R i^-t-zl 

O^. ( ^ O.U 31 i.O 


3.40 - 3.60 


16.3 ±2.5 ±1.9 


19.5 lb 5.1 ± 2.4 


3.60 - 3.80 


8.5 ±1.4 ±0.98 


8.1 lb 3.0 ± 0.96 


3.80 - 4.00 


3.67 ±0.25 ±0.42 


7.34 lb 0.77 lb 0.86 


4.00 - 4.20 


2.00 ± 0.13 ± 0.23 


3.89 lb 0.33 ± 0.45 


4.20 - 4.40 


1.242 ±0.078 ±0.14 


2.18 lb 0.19 =b 0.25 


4.40 - 4.60 


0.587 ± 0.043 ± 0.067 


1.20 lb 0.12 lb 0.14 


4.60 - 4.80 


0.409 ± 0.024 ± 0.047 


0.805 lb 0.073 ± 0.095 


4.80 - 5.00 


0.244 ± 0.017 ± 0.028 


0.502 ± 0.048 ± 0.059 




{j>b/{GeV/cY) 


\pu/\Kjt.v /O) ) 


5.00 - 5.25 


132.5 ±9.2 ±15 


Q41 _|_ 97 _|_ 40 


5.25 - 5.50 


70.2 ±5.6 ±8.2 




5.50 - 5.75 


36.3 ± 4.4 ± 4.3 


lU 1 31 IZ 31 ±0 


5.75 - 6.00 


17.2 ±2.2 ±2.0 


58.0 ±8.0 ±7.0 


6.00 - 6.50 


7.27 ±0.98 ±0.87 


26.3 ±3.6 ±3.2 


6.50 - 7.00 


2.63 ±0.43 ±0.32 


10.4 ±1.8 ±1.3 


7.00 - 8.00 


0.23 ±0.13 ±0.028 


3.02 ±0.64 ±0.38 


8.00 - 9.00 


0.035 ± 0.046 ± 0.005 


1.33 ±0.43 ±0.17 


9.00 - 10.00 




0.11 ±0.24 ±0.015 


10.00 - 12.00 




0.073 ± 0.050 ± 0.010 



TABLE IX. Invariant differential cross sections [Ed^a/dp^) for the inclusive reaction pp rjX. 
at 530 and 800 GeV/c, averaged over the rapidity intervals —0.75 < ycm < 0.75 and 



— 1.0 < T/cm < 0.5, respectively. 



PT 


pp @ 530 GeV/c 


pp @ 800 GeV/c 


(GeV/c) 


(nh/{GeV/cf) 


[nh/{GeV/c)'^) 


3.00 - 3.50 


22.6 ± 9.5 ± 2.7 


51 ±17 ±6.3 


3.50 - 4.00 


7.6 ±1.9 ±0.87 


11.4 ± 3.7 ± 1.4 


4.00 - 4.50 


1.19 ±0.14 ±0.13 


2.13 ±0.34 ±0.25 


4.50 - 5.00 


0.310 ± 0.035 ± 0.035 


0.678 ± 0.091 ± 0.080 




{pb/{GeV/cy') 


(p6/(Gey/c)^) 


5.00 - 5.50 


77 ± 12 ± 9.0 


220 ± 36 ± 26 


5.50 - 6.00 


23.7 ±5.0 ±2.8 


108 ± 16 ± 13 


6.00 - 7.00 


2.3 ±1.2 ±0.3 


20.6 ±4.8 ±2.5 


7.00 - 8.00 


0.41 ±0.31 ±0.05 


4.4 ± 1.8 ±0.6 


8.00 - 10.00 




1.05 ±0.70 ±0.14 
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TABLE X. The averaged invariant differential cross sections per nucleon [Ed^a / dp^) as a 



function of rapidity and px for the inclusive reaction pBe — > ryX at 530 GeV/c. The units are 
pb/(GeV/c)2. 



y 


cm 




PT (GeV/c) 










3.00 - 4.00 






4.50 - 5.00 


5.00 - 5.50 


-0.750 
-0.625 


- -0.625 

- -0.500 


1 Acinc\ _l_ nnnn _l_ i iic\n 

i4duu ± yuuu ± loUU 


1050 ± 310 ± 110 

yoU ± o2U ± iUU 




247 ± 62 ± 27 
318 ± 88 ± 34 


66 ± 25 ± 7 
96 ± 20 ± 10 


-0.500 


- -0.375 


24100 ± 4900 ± 2600 


1420 ± 280 ± 150 




374 ± 68 ± 40 


100 ± 26 ± 11 


-0.375 


- -0.250 


1730 ± 290 ± 190 




431 ± 60 ±47 


111± 18 ± 12 


-0.250 


- -0.125 


21200 ± 3800 ± 2300 


1930 ± 250 ±210 




486 ± 57 ± 52 


134 ± 19 ± 15 


-0.125 


- 0.000 


1750 ± 2U0 ± 190 




514 ± 56 ±56 


123 ± 17 ± 13 


0.000 


- 0.125 


26700 ± 3400 ± 2900 


1330 ± 170 ± 140 




482 ± 51 ±52 


121 ± 18 ± 13 


0.125 


- 0.250 


2U60 ± 160 ± 220 




364 ± 45 ± 39 


97± 15 ± 11 


0.250 
0.375 


- 0.375 

- 0.500 


22700 ± 3300 ± 2400 


1480 ± 120 ± 160 
139U ± 100 ± 150 




333 ± 39 ± 36 
327 ± 34 ± 35 


129 ± 20 ± 14 
94 ± 15 ± 10 


0.500 


- 0.625 


13200 ± 4100 ± 1400 


1220 ± 100 ± 130 




321 ± 36 ± 35 


90 ± 14 ± 10 


0.625 


- 0.750 


837 ± 74 ±89 




187 ± 28 ± 20 


56 ± 12 ± 6 






5.50 - 6.50 


6.50 - 8.00 








—0.750 
-0.625 


0.625 

- -0.500 


4.6 ± 7.0 ± 0.5 
13.4 ± 3.9 ± 1.5 


0.35 ± 0.33 ±0.04 








-0.500 
-0.375 


- -0.375 

- -0.250 


7.6 ± 6.9 ±0.9 
18.1 ± 4.0 ±2.0 


1.03 ± 0.46 ±0.12 








-0.250 
-0.125 


- -0.125 

- 0.000 


21.2 ± 4.7 ± 2.4 
19.4 ± 4.4 ±2.2 


1.59 ± 0.46 ±0.19 








0.000 
0.125 


- 0.125 

- 0.250 


21. 7± 4.3 ±2.4 
20.2 ± 4.0 ±2.3 


1.68 ± 0.45 ±0.20 








0.250 
0.375 


- 0.375 

- 0.500 


23.9 ± 3.8 ±2.7 
20.0 ± 3.9 ±2.2 


0.82 ± 0.41 ±0.10 








0.500 
0.()25 


0.625 
0.750 


18.6 ± .'^.3 ±2.1 
11.5 ± :-!.2 ± 


n.Gs = o.,s3 ± n.os 








TABLE XL 


The averaged invariant differential cross sections 


{Ed^a/dp^) 


as a function of 


rapidity and pT for the inclusive reaction pp 


r]X at 530 GeV/c. The units are 


pb/(GeV/c)2. 


Vcm 




PT (GeV/c) 










4.00 - 5.00 




5.00 - 6.00 


6.00 - 8.00 


-0.75 - 


■ -0.50 


520 ± 250 ± 60 




29 ± 19 ± 3 


1.6± 1.0 ±0.2 


-0.50 - 


-0.25 


700 ± 230 ± 80 




53± 16 ± 


6 


-0.25 - 


0.00 


1150 ± 190 ± 120 




57± 16± 


6 


1.1 ± 1.0 ±0.1 


0.00 - 


0.25 


740 ± 150 ± 80 




70 ± 18 ± 8 


0.25 - 


0.50 


800 ± 110 ±90 




46 ± 13 ± 5 


1.3± 1.2 ±0.1 


0.50 - 


■ 0.75 


573 ± 73 ±61 




48 ± 13 ± 5 
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TABLE XII. The averaged invariant differential cross sections per nucleon [Ed?a / dp^) as a 



function of rapidity and pt for the inclusive reaction pBe — > ryX at 800 GeV/c. The units are 
pb/(GeV/c)2. 



y 


cm 






PT (GcV/c) 












3.00 - 4.00 


4.00 


- 4.50 


4.50 - 5.00 


5.00 ■ 


- 5.50 


-1.00 
-0.875 


- -0.875 

- -0.750 


1 Qfinn _l_ 1 cnnn _l_ i nnn 
loUUU ± loUUU ± lUUU 


920 ± 
1360 ± 


580 ± 100 
620 ± 1.50 


— 

500 ± 240 ± 60 


116 ± 
87 ± 


47 ±13 
52 ± 10 


-0.750 
-0.625 


- -0.625 

- -0.500 


46000 ± 11000 ± 5000 


1980 ± 
1270 ± 


630 ± 220 
580 ± 140 


560 ± 130 ± 60 
540 ± 160 ± 60 


155 ± 

259 ± 


43 ±17 
48 ± 29 


-0.500 


- -0.375 


17300 ± 8100 ± 1900 


2330 ± 


550 ± 260 


640 ± 120 ± 70 


299 ± 


49 ± 34 


-0.375 


- -0.250 


3840 ± 


530 ± 420 


710 ± 150 ± 80 


253 ± 


56 ± 29 


-0.250 


- -0.125 


48200 ± 8800 ± 5300 


3520 ± 


570 ± 390 


570 ± 130 ± 60 


239 ± 


48 ± 27 


-0.125 


- 0.000 


4130 ± 


520 ± 460 


940 ± 140 ± 100 


292 ± 


53 ± 33 


0.000 


- 0.125 


34500 ± 6700 ± 3800 


4070 ± 


600 ± 450 


1200 ± 160 ± 130 


311 ± 


64 ±35 


0.125 


- 0.250 


2840 ± 


510 ± 310 


1220 ± 180 ± 140 


386 ± 


76 ±44 


0.250 


- 0.375 


29000 ± 7100 ± 3200 


3670 ± 


350 ± 410 


1170 ± 130 ± 130 


275 ± 


57 ±31 


0.375 


- 0.500 


2230 ± 


260 ± 250 


1090 ± 110 ± 120 


282 ± 


48 ±32 






5.50 - 6.50 


6.50 


- 8.00 


8.00 - 10.00 






-1.00 
-0.875 


- -0.875 

- -0.750 


9 ± 10 ± 1 
22.5 ± 8.9 ±2.6 






1.50 ± 0.93 ±0.19 






-0.750 
-0.625 


- -0.625 

- -0.500 


51± 16±6 
41± 11±5 


6.8 ± 


2.0 ±0.8 


0.52 ± 0.93 ±0.07 






-0.500 
-0.375 


- -0.375 

- -0.250 


34 ± 12 ± 4 

51± 11 ±6 


4.3 ± 


1.6 ±0.5 


0.55 ± 0.52 ±0.07 






-0.250 
-0.125 


- -0.125 

- 0.000 


86± 16 ± 10 
79 ± 13 ± 9 


6.4 ± 


1.5 ±0.8 


0.68 ± 0.23 ±0.09 






0.000 
0.125 


- 0.125 

- 0.250 


93 ± 17± 11 
70 ± 18 ± 8 


8.3 ± 


2.2 ±1.0 


0.61 ± 0.28 ±0.08 






0.2-50 

0.:-!75 


0.375 
0.500 


61 ± 15 ± 7 

51= i:-! ± (i 


7. 1± 


l.S ± 0.9 


0.1S± 0.27 = 0.00 







TABLE XIII. The averaged invariant differential cross sections clS cl function of 



rapidity and pT for the inclusive reaction pp rjX at 800 GeV/c. The units are pb/(GeV/c)^. 



ycm 




PT (GeV/c) 






4.00 - 5.00 


5.00 - 6.00 


6.00 - 8.00 


-1.0 - -0.75 


1060 ± 430 ± 120 


73 ± 50 ± 8 


6.9 ± 3.6 ±0.8 


-0.75 - -0.50 


360 ± 520 ± 40 


97 ± 49 ± 11 


4.9 ± 7.1 ±0.6 


-0.50 - -0.25 


1930 ± 420 ± 210 


194 ± 40 ± 22 


14.4± 7.6 ± 1.7 


-0.25 - 0.00 


1880 ± 440 ± 210 


256 ± 47 ± 29 


22.3 ± 6.0 ±2.6 


0.00 - 0.25 


1980 ± 470 ± 220 


167 ± 55 ± 19 


16.9 ± 7.6 ±2.0 


0.25 - 0.50 


1230 ± 250 ± 140 


198 ± 50 ± 22 


9.6 ± 5.0 ±1.1 
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